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NOTATIONS 


x,v 

z 

i 


R,S 


L 

K 

n,m 


XX,  yy 


xy 

U 


U_ 


7^ 


Hectangular  coordinates  in  the  z plane. 
Complex  variable  z = x + iy. 


i = 


(-1) 


1/2 


Polar  coordinates  in  the  w plane. 

Complex  variable  w = Re^®. 

Boundary  of  the  hole  in  the  z plane. 

Unit  circle  in  the  w plane. 

Direction  cosines  of  the  outer  normal. 
Components  of  surface  forces  per  unit  area 
at  a point  of  the  boundary. 

Normal  components  of  stress  parallel  to 
the  X-  and  y-axes. 

Shearing  stress. 

Biharmonic  stress  function  U = U(x,y). 
Subscripts  x and  y denote  partial  differ- 
entiation with  respect  to  the  indicated 
variable . 


\7^  = ff 


3x^  3y^ 


(i) 


CHAPTER  I 


INTRODUCTION 

This  is  a mathematical  analysis  of  the  stress  distri- 
butions existing  in  a large,  uniformly  loaded  isotropic 
plate  weakened  by  a small  hole.  We  will  consider  the 
plate,  whose  plane  is  the  xy_nlane,  to  be  subject  only  to 
forces  directed  parallel  to  the  x-  and  y_axes.  It  has 
been  shown  that  the  solution  to  such  problems,  when  body 
forces  are  absent  or  constant,  reduces  to  the  integration 
of  the  differential  equation^ 

U = 0,  (1) 

having  regard  to  the  boundary  conditions 

- 1 
= nxx  + mxy  1 

V - - r 

The  stress  comnonents  are  represented  in  terms  of 
the  biharmonic  stress  function  by  the  set  of  equations 

Sc  = Uyy.  fy  = iy  = -U^,  (3) 

where  subscripts  on  U denote  partial  derivatives  with 
respect  to  the  indicated  variable. 

The  boundary  curve  L of  the  hole  in  the  z plane, 
will  be  taken  to  be  given  in  terms  of  the  nolar  angle  Q 

§ See  reference  5 of  Bibliography. 
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in  the  w plane,  by 


X = r ' 

Cos 

5S 

+ 

r" 

Sin 

53 

+ 

t ' 

Cos 

23 

+ t" 

Sin 

2S 

+ 

8 ' 

Cos 

3 

- 

s" 

Sin 

3, 

II 

Cos 

53 

- 

r ' 

Sin 

53 

+ 

t" 

Cos 

23 

_ t' 

Sin 

23 

+ 

s" 

Cos 

3 

+ 

s ' 

Sin 

3. 

V (4) 


The  granh  of  these  equations  for  various  values  of  r» , 
r",  t',  t”.  S',  and  s",  will  be  discussed  later. 

We  use  the  mapping  function 

z = ■p(w)  » rw~  + tw“  + sw,  (5) 

to  man  the  region  exterior  to  L onto  the  region 
exterior  to  the  unit  circle  F in  the  w plane,  where 

r = r*  +ir",  t = t'  + it",  s = s'  + is".  (6) 


To  find  the  biharmonic  stress  function  U we 
formulate  the  boundary-value  problem  in  the  form  of 
integral  ecuations^,  in  order  to  anply  the  method  most 
often  associated  with  the  name  of  N.  I.  MusheliSvili . 

^e  has  shown  that  if  P(w)  is  rational,  then  the 
solution  of  the  integral  equations  reduces  to  evaluating 
integrals  of  Cauchy's  type. 


5 See  reference  4 of  Bibliography. 

See  references  7 and  8 of  Bibliography 


CHAPTER  II 


THE  HOLE 

In  our  discussion  of  the  hole  whose  boundary  L in 
the  comnlex  z plane  is  given  by  Eqs.  (4),  we  first  con- 
sider the  conformal  mapping  of  the  complex  w plane  unon 
the  z plane.  For  the  mapping  to  be  conformal  a theorem 
which  we  shall  call  Theorem  I must  hold. 

Theorem  T.  At  each  point  where  a function  P(w)  is 
analytic  and  P*(w)  ^ 0,  the  mapping  z = P(w)  is  con- 
formal. 

In  order  that  the  unit  circle  in  the  w plane  (here- 
after called  P)  shall  map  upon  L,  let 

z = P(w)  = rw~^  + tw“^  + sw.  (7) 

In  the  apnlications  which  follow,  attention  is  confined 
to  those  problems  for  which  the  singularities  of  P(w) 
and  the  zeros  of  P*(w)  lie  on  the  boundary  or  outside 

of  it.  Prom  Eq . (7)  we  see  that  the  point  w ^ oo 

corresponds  to  the  point  at  infinity  in  the  z plane, 
and  that  we  have  no  poles  at  points  in  the  date.  We 
note  that  P(w)  is  a single-valued  analytic  function 
everywhere  in  the  plate.  The  inverse  function  w = g(z) 
is  a multiple-valued  analytic  function,  and  for  any 
value  of  z there  are  in  general  one  or  more  values  of 
w.  Hence  we  find  it  necessary  to  assign  such  values  to 
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r,  t,  and  s,  that  the  branch  points  in  the  z plane 
shall  lie  on  the  boundary  or  outside  of  it,  since  in 
the  neighborhood  of  a branch  point  the  manping  by 
means  of  an  analytic  function  ceases  to  be  conformal. 

A criterion  for  locating  such  points  is  given  by  the 
following  theorem. 

Theorem  TI.  Given  a multiple-valued  analytic 

function  w = g(z),  whose  inverse  function  z = T’’(w)  is 
single-valued.  If  the  derived  function  ?'(w)  has  a 
zero  of  order  h,  at  a finite  point  w^,  or  a nole  of 
order  h + 2,  then  g(z)  has  a branch  point  of  order  h 
at  the  corresponding  point  z^.- 

v/e  can  also  state  that  if  i^*(w)  has  a pole  of  order 

h at  w = 00  , then  g(z)  has  a branch  point  of  order  h at 

the  corresnonding  point  in  the  z plane. 

Tn  annlying  Theorem  TI,  we  do  not  need  to  consider 
the  poles  for  w = 0,  since  we  will  alwavs  use  values  of 
w so  that  |w|  =1,  in  order  to  stay  in  the  proner 

branch  of  the  w function.  ?rom  Ka . (7) 

■p*(w)  = - 5rw"^  - 2tw~^  + s.  (8) 

To  determine  the  zero  points,  we  have 

§ ■por  proof  of  Theorem  II,  see  page  335  of  reference 
6 of  Bibliography. 
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where 


'’’he  six  zero  points  are 


and 


k = (2t  + + 20rs 


(9) 


where  1,  p,  and  p^,  are  the  cube  roots  of  one. 


^1  = k,  = pk,  w^ 


(10) 


The 


branch  points  are  obtained  bv  substituting  the  values 
of  w froic  Eqs.(lO)  into  Eq.(7).  Parameters  r,  t,  and 
3,  will  be  taken  so  that 


’^Hien  |j(  - 1,  or  Ikl  =1,  the  boundary  curve  will  have 
cusps,  as  in  Pig.  2,  Appendix  T.  Since  these  branch 
points  are  the  zeros  of  P*(w),  inspection  of  Eqs.fSg) 
to  will  show  that  the  stresses  become  infinite  at 

these  points.  In  the  neighborhood  of  such  points  a 
plastic  flow  win  occur  which  will  eliminate  the 
infinite  stresses,  and  our  equations  are  not  valid  in 
such  a neighborhood.  When  |j|  and  |kl  <1,  the 

boundary  curve  has  continuous  curvature  and  our  solution 
is  valid.  The  mapping  z = F(w)  is  then  conformal  for 
all  values  of  |w|  = 1. 


IJI  = 1,  Ik|  ^ 1. 


(11) 
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We  now  consider  some  illustrations  of  various 
holes.  The  simplest  hole  is  for  r=t=0,  s=s*, 
which  is  a circle.  This  problem  has  been  previously 
solved  so  we  omit  this  case.^ 

CASE  ONE.  let  r = 0,  and  let  t and  s be  real. 

Eouations  (4)  reduces  to 

X = t Cos  2B  + s Cos  S, 
y = - t Sin  2B  + s Sin  S. 

Erom  Eq.  (8) 

E* (w)  = _ 2tw“^  + 3.  (13) 

'^he  zero  points  of  E*(w)  = 0,  are 

Wi  = (2ts"^)^'^^,  w^  = pw^,  w^  = p^w^.  (14) 

Erom  Eq . (7) 

z = p(w)  = tw“^  + SW.  (1‘i) 

On  substituting  the  values  of  w from  (14)  into  (15), 
we  obtain  the  branch  points 

= |(2ts^)^/5,  = pz^,  z^  = p^z^.  (16) 

In  calculating  the  above  values  of  z^,  z^,  and  z^,  we 
have  made  use  of  the  relations 

^ (p2)3  ^ ^ 

which  follow  from 


§ See  reference  5 of  Bibliography. 
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p . 1 * p2  - 1 

2 2 ,p--2~2 

Setting  t = 1,  s = 5,  in  Eq.(l2) , we  have 

X = Cos  2S  + 5 Cos  0, 

V = - Sin  20+5  Sin  0. 

'T’he  curve  obtained  from  Eqs.fl?)  is  given  in  Appendix  I, 
■Pig.  1.  The  branch  points  calculated  from  Eqs.(16)  give 
7,^  = 5.526,  Zp  = - 2.763  + 4.786i,  = _ 2.763  -4.786i, 

which  lie  inside  L. 

Erom  Eos. (12),  on  setting  s =2,  t =1,  we  have 
X = Cos  20+2  Cos  0, 
y = - Sin  20+2  Sin  0. 

■'J'rom  Eq . (8) 

?*(w)  = _ 2w“^  + 2.  (19) 

T’he  zero  points  of  E*(w)  = 0,  are 

w^  =1,  = p,  w^  = p^.  (20) 

Substituting  the  values  of  w from  Eqs.(20)  into 
Eos. (16)  we  obtain  the  branch  points 

^1  =5,  = 3p,  z^  = 3p^.  (21) 

The  branch  points  in  Eq.(21)  lie  on  I,  since  k = 1 in 

Eq.(ll).  The  curve  obtained  from  Eqs.(l8)  is  given  in 
Appendix  I,  Pig.  2. 
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CASE  TWO.  Let  t = 0,  and  let  r and  s be  real. 
Eos.  (4)  reduces  to 


X = r Cos  5S  + s Cos  6, 


y = - r Sin  53  + s Sin  3. 


(22) 


Erom  Ea . (8) 

E'(w)  = _ 5rw"^  + s«  (23) 

The  six  zero  points  of  E*(w)  ■ 0,  are  obtained  from 
Eos. (9),  where 

k * (5rs"^)^'^^,  j • « k. 


From  Bq. (7) 


z = E(w) 


+ sw. 


(24) 


On  substituting  these  six  values  of  w into  Eq.(24), 

we  obtain  the  branch  points 

2 ? 
z^=  a,  82=  pa,  z^=  p a,  z^a  _ a,  z^c  _pa,  Zg«  «p  a,  (25) 


where 


.576 


(26) 


Setting  8 = 20,  r = 1,  Eqs.(22)  becomes 


X 


y 


Cos  53  20  Cos  Q, 

Sin  53  + 20  Sin  80 


(27) 


The  curve  obtained  from  Eqs.(27)  is  given  in  Appendix  I, 
^i^:.  3.  We  note  that  the  curve  is  in  the  shape  of  a 


« 
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hexaf^on  where  the  corners  have  been  rounded  off. 
Substituting  s = 20,  r = 1,  into  Eq.(26),  gives  a = 
(19.044).  Bv  inspecting  Eos. (25),  we  see  that  the 
bra’^ch  points  lie  on  a circle  of  radius  a,  at  intervals 
of  60°,  and  are  inside  of  L. 

Let  r = ],  s = 10.  From  Eqs.(22),  we  obtain 


X = Cos  5S  + 10  Cos  8, 


y = - Sin  58  + 10  Sin  8. 


(28) 


The  curve  obtained  from  Eqs.(28)  is  given  in  Appendix 
I,  Fig.  4.  All  the  branch  points  lie  inside  I.  There 
are  no  sharp  corners  on  the  curve. 


CASE  TT^^EE.  This  case  includes  all  tv^e  holes 
obtainrble  from  Eos. (4),  which  are  not  already  included 
in  Case  T or  Case  IT . The  stress  components  given  in 
Eos. (93)  to  (95),  and  the  biharmonic  stress  function  TJ, 
Eo . ( 87 ) , are  derived  for  real  or  complex  values  of  r, 
s,  and  t.  We  shall  give  two  examnles  of  holes  that 
can  be  obtained  from  Eqs.(4).  Let  r*  = r"  =1,  t»  = 
t"  = 0,  s'  = s"  = 10.  Then  Eqs.(4)  reduces  to 

X = Cos  58  + Sin  33  + 10  ("-os  8 - Sin  8),] 

, I 

y = Cos  S8  _ Sin  58  + 10  (Cos  0 + Sin  8).J 

'^he  hole  whose  boundary  is  given  bv  Eos.  (29''  is  given 
in  Appendix  I,  ‘'^'ig.  , Equation  (7)  reduces  to 
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z = F(w)  = (1  + i)w”^  + 10(1  + i)w.  (30) 

Hence 

P*(w)  = - 5(1  + i)w~^  + 10(1  + i),  (31) 

and  It  is  seen  froir  Eq.(31),  that  the  zero  noints  of 

F»(w)  = 0,  are  = 1/2.  Thev  lie  inside  and  the 
branch  ooints  lie  inside  L. 

Let  r*  =r”  =1,  t*  =t”  =0,  s*  =s"  =20. 
Ecuation  (4)  reduces  to 

X = Cos  5B  + Sin  5S  + 20(Cos  Q _ Sin  S), 

y = Cos  50  - Sin  56  + 20(Cos  0 + Sin  0). 

In  Appendix  I,  Pig.  6,  we  have  the  curve  obtained 
from  Eqs.(32).  It  is  similar  to  the  curve  in  Pig.  3. 


CHAPTER  III 


the  stress  PIR^GTION 

The  general  blharmonlc  stress  function  can  be 
expressed  in  ter«s  of  t«o  analytic  functions  of  the 
co.ulex  variable  z.  let  . p. 

7‘n  = v^V^u  = V^P  = 0. 

hence  P is  a harmonic  function.  The  function  P being 
harmonic,  „e  can  construct  an  analytic  function 

. Pfx.y)  + lQ(x,y), 

where  the  functions  P and  0 sati  q-Ptr  +v,  n 

^ satisfy  the  Cauchy.Riemann 

eouations 

~ y » P = _ Q , 

let  us  define  an  analytic  fuLtlon  ^(s)  by  the  formula 
then  = P . Iq  = 1/4 ^(z)dz, 

0'(^)  = . 1,^  = i/4p(,)  ^ ^ 

and  since  0=0  n 

^ ~ ~ that 

= - 1/4Q. 

It  then  follows  that 


- xp  _ yq)  = 0, 

- that  the  hihar.onic  function  H has  the  for. 


U = xp  + vq  + p^, 


where  p^  ig  go^jg  harmonic  function. 
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Let  us  construct 

•z 


the  function 

z)dz  = + iq 

0 


o’ 


where  is  the  harinonic  conjugate  of  p^. 
function  U,  in  terms  of  j?5(z)  and  T(z),  is 

IT  = Re(z0(z)  + I T(z)dz) , 


The  hiharmonic 
then  given  by 

(33) 


0 

where  z = x _ iy,  and  the  symbol  Re  denotes  the  real  part 
of  the  expression  following  it.  As  a consequence  of  the 
last  equation,  the  stresses  can  be  exnressed  in  terms  of 
0(z)  and  T(z).  The  calculations  of  the  stress  components 
give  the  following  equations. 


= Re(_  i0(z)  + iz0'(z)  4 iT(z)).  (34) 

^ = '^e(20*(z)  _ zj?f"(z)  _ T'(z))  = 

= /’(z)  4 0*(z)  _ ■|(z/"(z)  4 z^"(z)  4 T’(z)  4 T'(z)). 

= Re(i?^(z)  4 z0*  (z)  4 T(z) ) . (36) 

fy  = Re(2j?^' (z)  4 z/"(z)  4 T'(z))  = U ^ 

XX 

= 0^{z)  4 0'(z)+  ~(z0"(z)  4 Z0"(z)  4 T'(z)  4 T'(z'>'). 

XV  = - = _ Re(izj?^"(z)  4 iT'(z'>) 

= - il(z0"(z)  4 T*(z)) 

- - i-|-(z/"(z)  _ z/"(z)  4 T*(z)  _ T»(z)). 


(38) 
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■ny  irs-oecting  ■^os.(3‘^)  and  (37),  we  see 

Q + yy  = 4Pe(i?^'(z))  = 20*{z)  + 20'(z).  (39) 

■^rom  Eqs.(35),  (37),  and  (38),  we  obtain 

-■Q  + ^ + 2iQ  = 2z/"(z)  + 2T*(z).  (40) 


Consider  an  arc  A?  of  some 
curve  in  the  nlane  of  the  plate, 
e.o;  the  boundary  curve, 

We  shall  orient  the  positive 
direction  of  the  norm.al  F and 
of  the  tangent  T as  indicated  in 
■^ig.(a).  Let  X^t<3s  and  X^^ds  represent 
the  X-  and  y_components  of  the  force 
acting  or  the  element  of  arc  ds  of  AE, 
from  t^e  positive  side  of  the  normal  N.  'Ve  have  from 
Eqs. (2) 


= XX  Cos(x,’^0  + ^ Gos(y,>’) 

= Tj^yCos(x,^0  - n^^Cos(y,^0, 

^ Cos(x,N)  + yy  Cos(y,M) 

= n^/os(y,N)  _ U^yCos(x,N). 

Eut 

Cos(x,N)  = Cos(y,T)  = dy/ds, 
Cos(y,Y)  = _ Cos(x,T)  = - dx/ds, 
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so  that 

= d/ds(lT^),  = - d/ds(U^).  (43) 

■'^’orrain^  the  expression 

+ iY^  = d/ds(Uy  - iU^)  = - id/ds(U^  + iUy), 

we  obtain 


(X^,  + iY^)ds  = > id(U^  + iUy). 
Then  the  quantity 


iY^)ds  = 


f" 

ijd(U^  + iUy)ds 


+ z0* (z) 


(44) 


- -n® 

+ T(z)  , 
— A 


represents  the  resultant  force  on  the  arc  AB,  where  we 
have  substituted  the  values  of  IT^  and  from  Eqs.(34) 
and  (3^).  Equation  (44)  expresses  the  boundary  con- 
ditions in  terms  of  0(z)  and  T(z),  when  forces  are 
prescribed  on  a boundary  curve.  In  the  case  at  hand, 
our  plate  is  subject  onlv  to  stresses  at  infinity.  Let 
the  prescribed  forces  be  (xx)_=  S,  , (yy)  = S«,  (^) 

= ^22*  ^®^ce,  no  forces  act  on  the  boundary  of  the 

hole,  and  the  left  hand  side  of  Eq.(44)  becomes  zero. 

We  shall  find  it  more  convenient  to  deal  with  these 
boundary  conditions  in  the  w plane.  Becalling  Eq.(‘^) 

z = P(w)  = sw  + tw~^  + rw“^,  (45) 
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we  have 


z*  = F*(w)  = s - 5rw"  . 


(46) 


For  the  sake  of  simplicity,  we  introduce  the  following 
notations ; 


Then 


and 


j2^(z)  - iZl(F(w) ) - 0(w) , 
T(z)  - T(F(w) ) - T(w) . 


0*{z)  = (w)/F* (w) , 


T * ( z ) = T ' (w)/F  * (w) . 


The  value  of  w on  K shall  he  denoted  hy  v.  We  have 


iS 


_ -1 


V = e-"“  = (v)  . 

Noting  that  F corresponds  to  L,  on  the  boundary 
z0' (z)  = F(v)0* (v) 

, 0.,Z) 

s _ 2tv”^  _ 5rv“" 


V (sv  + tv  -t-  r) 
s _ 2tv^  -5rv^ 


0*  (v) 


(47) 

(48) 

(49) 

I 

(60) 


(51) 


= H(v)0» (v) . 

Substituting  Fq.C^l)  into  Eq.(44),  the  boundary 
conditions  read 


0(v)  + H(v)0*(v)  + T(v)  = 0, 


(62) 
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or  - - _ « 

0(y)  + H(v)j?f»(v)  + T(v)  = 0,  (53) 

where  Eq.(53)  is  the  conjugate  of  Eq.(52),  and 
H(v)  = v^(s  + tv^  + rv^)(sv^  _ 2tv^  - 5r)"^. 

It  is  known  that  when  there  is  no  unbalanced 
force  acting  on  the  hole  and  when  the  stresses  at 
infinity  are  uniform,  the  functions  0(w)  and  T(w) 
take  the  form^ 

00 

J^(w)  = Cw  + ^ » (54) 

n=l  " 

00 

T(w)  = (D  + iD*)w  + 21  ^ (55) 

nfl 

where  the  constants  C,  T),  and  are  real,  and  the 
constants  a^^  and  b^  may  be  real  or  complex.  To  find 
the  value  of  the  constant  C,  we  aoolv  Eqs.(39),  (45), 
(48),  and  (54). 

(^)qo+  (^)qo  = (^  + ^)qo  = 2(j2(»(2:)  + R’(z))^ 

= 2(i?f' (w)/E*  (w)  + (w)/e*  (w)  )^ 


00 

= 4Cs*(ss)“^  = + S^. 


= 2 


00  - 
C - ^ na^w"’^“ 

s j-2tw”^  -5rw~^ 


c -r 


C30  _ _ 


n=l 


nanW 


-n-1 


s -.2tw”^  -5rw”^ 


§ See  reference  3 of  Bibliography. 


17 


To  calculate  the  values  of  7}  and  D*  in  Eq.(5‘=;) 
have  froir  Eq  . ( 40) 


^^^00  - ^^^00  + 2iS 


12 


= 2(z^"(z)  + T» (z))^  . 


Ta>lr.g  the  derivative  of  Eq.(48),  we  obtain 

0"(z)dz/dw  = fE*  (w))"^[j*  (w)>^"(w)  - ;?*(w)E"(w) 


Eence 


0"(z)  = (E*  (w)  )”^  jp» 'w)jZ^"(w)  _ /*(w)E"(w^. 

Substituting  the  values  of  (w) , 0"(w),  E'(w),  and 
E’'(w),  we  have 

•T  c 00  ^ 

(s  - 2tw“  _ 5rw“  )(  n(n  + l)a^w~^”^ 
^"(z)  


(s  _ 2tw”^  _ 5rw"^)^ 

00  ^ A 

(C  - ^ na^w”’^"-'-)  (6tw"^  + 30rw“') 


(s  - 2tw”^  - Srw**^)^ 


Bv  inspecting  Eq.(S8),  we  see  (/"(z))^  = 0.  Fence 
Ec.(56)  reduces  to 


2(T» (z) ) 


00 


00  - 

2(D  + iD»  - ^ nb^w-"-^) 


s - 2tw~^  _ Srw”^ 


00 


.-1, 


= 2s”  (D  + iD»)  = -3^  + 

Let  B = s(S^  + S^)(4s*)“^,  and  B’  = 1/2(8^  - 3^). 


, we 

(56) 

]• 

(57) 
’) 

(58) 
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Then  Eq3.(54)  and  become 

00 

= sEw  + ^ a^w“  , (59) 

°° 

T(w)  = s(B»  + 

The  actual  process  of  solving  the  problem  consists 
of  multiplying  Eqs.(52)  and  (53)  by  (2TTi(v  _ w))‘'^dv 
and  integrating  around  IT.  The  inteprration  will  be 
performed  with  the  aid  of  the  following  two  theorems. 


Theorem  T!^T . 


If  f(w)  is  an  analytic  function  within  K, 


except  perhaps  at  w = 0,  where  it  has  a pole  with  a 

N 

nS 


principa]  part  of  then  (^f  (v)  ( 2tt  i(v-w)  )“^dv 


n=. 


A^w”’^,  for  Iwl  > 1 . 


To  prove  T'heorem  III,  we  represent  f(w)  by  its 
Laurent  series 


00 

y B w 

— m 


m=o 


m 


in  a neighborhood  of  the  nole  w = 0.  Then  applying 
the  Residue  Theorem,  we  have 


~^dv  =-(2TTi)“ 


z;  vPw-p-idv  = 

p=o 
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Theorem  IV. 

If  f(w)  is  an  analytic  function  outside 
of  IT,  except  perhaps  at  w = oo  , where  it  has  a pole 
with  a principal  part  of  a P +v, 

A W f Xii©n 

p»o  ^ 


f 


f(v)(2TTi(v  - w))“^dv  = _ f(w) 


ZApwP, 
p=o  ^ 


for  w >1. 


To  prove  Theorem  IV,  we  represent  f(w)  by  its 
Laurent  series 


N 


.P  + 


p=o 


00 


m=. 


m 


.~m 


in  a neighborhood  of  the  pole  w = oo . Then,  applying 
the  Residue  Theorem,  we  have 


^f(v)(2TTi(v  - w))-^dv  = - (27Ti)-^^(v)  f^v^-^-^dv 


N 

= - f(w)  J2.  A wP,  for  )wl>l. 
p=o  P 
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let  us  first  multiply  Eq.('52)  by  ( 2TT i( v-w)  )”^dv 
and  integrate  around  K.  We  have 


dv  ^H(v)?t(v)dv  ^T(v)  dv 
J2TTi(v-w)  J 2TTi(v-w)  J^2TTi(v_w 


) 


= 0. 


Applying  Theorem  III  to  the  first  integral  in  the  last 
equation,  we  have 


(2TTi) 


-1 


00 


n=. 


a v-"") 
n ' 


dv 


V _ w 


00 


(61) 


= - ^ ^n’^”  ~ - jZ^(w)  + sBw. 


n=l 


To  evaluate  the  second  integral,  we  shall  write  H(v) 
in  the  following  form. 


H(v)  = r(s)“^v”'’  + Ev"^  + Hq(v),  (62) 

where  the  exponents  of  v in  H^(v)  are  all  positive, 

and  _ _ _ 

E = (ts  + 2tr)(s)“^. 


We  note  from  Theorem  IV,  that  the  value  of  the  Integral 
contains  only  negative  powers  of  w.  Hence 


(2TTi)“^ 


P(v)^*(v) 

V - w 


dv 


= (2TTi)~^ 


+ Ev“^)(sB 


^na  v’^'*^^)  (v_w)”^dv 
n=l  " 


BEsw“^ 


+ 


na^w 

n 


(63) 
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Applying  Theorem  III  to  the  third  integral,  we  have 


"T(  v)dv 


= (2TTi)"^  J^(v-wf(sB*v”^  - iaS^2''^~^  + ^ ^n^” 


= - s(B«  - 


(64) 


Adding  the  results  of  Eqs.(61),  (63),  and  (64),  since 
their  sum  is  zero,  we  have 

-0(w)  + sBw  _ Brw”^  « BsEw""^  + r(s)"^ 

- s(B'  - = 0.  (65) 


? na^w 
nfr 


Equation  (65)  determines  0(w)  except  for  the  constants 
a^,  a2,  a^.  These  coefficients  are  found  by  multiply- 
ing  Eq.(52)  successively  by  (2TTi(v  _ w'l)”^  v dv, 
(27Ti(v  - w))“^  v^  dv,  and  (2TTi(v  - w))~^  v^  dv, 
and  repeating  the  integration  around  K.  We  have 


f^0(v)dv  . i!vH(v)0»(v)dv  ^ ^vT(v)  dv 

2TTi(v-w)  ^ ^ 2rri(vlw)  ^^Tri(v-w) 

Inspection  of  Eq.(64)  shows 


0.  (66) 


0. 


(67) 
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InsTiectlor  of  ?lq.(6l)  shows 


^ “V 


Inspection  of  Eq.(63)  shows 


Upon  substitution  of  these  values,  Eq.(66)  reduces  to 


-wj^(w)  + sBw^  + a,  -Erw-^  -BsEw-^  4 — J^na  w’^~^=0.  (68) 

; n.l  " 

Multiplying  Eq.(''2)  by  (2TTi(v  - w))“^  v^  dv,  and 
applying  the  results  of  Eq.(67),  we  have 


^x5/(”l4w  + . 0 

^’tTl(v_^  ^ ^ 2^i(vlw)  * 

The  value  of  the  first  integral  is  easily  seen 
from  Eo.(61).  We  have 


(69) 


/( v)dv  = - a^w^^"  = _w‘^0(w)  4 sBw^  4 a, w 4 a«. 


27T  i ( v_w)  n=3 


The  second  integral  is  easily  evaluated  with  the  aid 
of  Eo . ( 63 ) . We  have 


(27Ti(v-.w'' )”^v^H(v)i?J' (v)dv  = _Rrw“^  4 r(s)”^a^w“^. 


,-3 


=r^-l^  _-l 
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Fence,  'F]q.(69)  reduces  to 


-W  ^(w)  + sBw^  + a^w  + -Brw~^  + r(s )”^a^w”^=0. 

Multiplying  Eq.(52)  by  (27Ti(v  - w))”^  dv, 
applying  the  results  of  Eq.(67),  we  have 


00 


i' 

f 

f 


l(27Ti(v  - w))-^  v^  0{v)  dv 
'2TTi(v  _ w))"^  v^hCv)?*  (v)dv  = 0. 
Inspecting  Eq,(61),  we  see 

'2TTi(v  _ w))“^  v^  0{v)  dv 

UU  <2  0m 

= - & _ w^j2<(w)  + sEw^  + a,w^  + a^w  + 

n=4  12 


Inspecting  Eq.(63),  we  see 
^(2TT±(v  - w) )“^v^H(v)^» (v)  dv  = . Erw”^. 

Hence  Eq,(71)  becomes 


- w^0(w)  4 sEw'^  4 a^w^  4 a^w  4 a^  _ Erw”^  = 0, 


(70) 

and 


(71) 


(72  ) 


Let  us  rewrite  Eqs.(65),  (68),  (70),  and  (72). 
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ra^(s)"’^w"^  + 2ra2(s)“^w”^  + 3ra^( 

= ^^(w)  _ sBw  + rBw"*^  + BsBw“^  + 3(7?’  _ 

“**T1 

a^  + ra^(s)“  w“  + 2ra2(s)“  w""^ 

m w0(w)  - sBw^  + rBw*"^  + BsKw"^.  (74) 

a^w  4 ra^(3)”^w”^  4 

= w^^(w)  _ sBw^  4 rBw**^.  (73) 

2 

■+■  gl  + ^3 

= w^0(w)  - sBw^  4 rBw~^.  (76) 

Multiplying  Eq,(73)  by  w and  subtracting  from  Eq.(74), 
we  obtain 

a^  - 3ra^(s)”^  = - s(B»  - iS^^) . (77) 

Multiplying  Eq.(75)  by  w and  subtracting  from  Eq,(76), 
we  obtain 

a^  = ra^(s)”^,  or  a^  = ra^s"^.  (78) 

Substituting  the  value  of  a^  into  Eq.(77),  we  have 
= s(s)^(B‘  - iS^2)^5rr  . ss)"^. 


(79) 
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From  Eqs.(78)  and  (79)  we  obtain  the  value  of  a^. 


= rss(B*  + 13^2^  * 


"'-1  (80) 


Multiplying  Eq.(74)  by  w and  subtracting  from  Eq.(75), 


a^  - 2ra2(8)“^  = - EsE, 


which  gives 
^2 


a«  = (4rr  - ss)”^  (2rs^BE  + s(s)^BE),  (81) 


and  0(w)  is  completely  determined.  Setting  a^  c 0, 

and  ae  = - rB,  from  Eq.(76)  we  have 

0 


0(w)  = sBw  + 


ni 


_-n 


(82) 


To  determine  T(w)  we  multiply  Eq.(53)  by 


(2‘rTi(v  - w))”  dv,  and  integrate  around  K. 

i?i(v)  dv  + i?H(v)i2f' (v)dv  + ^ T(v)  dv 
^?Vl(v-W)  ^ 2^tv-^  f 2-rfi(v_w 


) 


= 0.  (83) 


Performing  the  first  integration  with  the  aid  of 
Theorem  III,  we  obtain 


o P(sBw“^  + ^ • 

- j>  2mUi  w)  ^ 


dv  = - sBw”^, 


Performing  the  second  integration  with  the  aid  of 
Theorem  IV,  we  obtain 
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H(w)^*(w)  + 21 

p=o  P 


where  we  need  to  calculate  the  values  of  Ap,  which  are 
the  coefficients  of  positive  powers  of  w in  the 
expansion  of  H(w)j2l»(w),  We  have 


H(w)i?f*(w)  = 


(rw^^  + tw^  + sw^) 


(sw^  _ 2tw^  - 5r) 


( sB 


> £ 
n=l 


na^w  ) 


= (rs-V  + Ew^  + Hq(w))(sB  - 2^  na^jW-"-^), 

where  H^(w)  consists  of  terms  containing  w with  nega- 
tive exponents.  Noting  that  we  may  drop  Hp(w)  and 
a^w”  , we  have 

(r8"^w^  + Ew^)(sB  - a^w”^  _ 2a2w“^  - 3a.^w“^) 

= rBw'  _ rs’^a^w^  + (BsE  _ 2rs~^a2)w^ 

- 3ra^s"^w  - Ea^  - E(2a2w“^  + 3a^w"^). 


The  values  of  Ap  are  then  obtained  from  the  last 
expression. 


27 


^2 


^3 

% 

B 


B* 


E 


2 - - P-  - - T >(84) 

= (2rs  BE  + s(s)'^BE)(4rr  - 3s)“'^. 


rss(B»  + - sa)“^. 


= 0. 


= - rB. 


s(Sj  + S2)(4s’)“^ 


(S^  - S^)/2. 


/ 


i 

/ 


= (s)“^(ts  + 2tr). 


Applying  Theorem  ITT  to  the  third  integral  in  Eq . (83), 
we  have 


§■ 


T(v)  dv  : 
27Ti  ( v_w) 


/ 


00 

(sB*v  + iaS^^"''’  + ^ 

(2TTi(v  - w) 


b„v-") 


dv 
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00 


n=j 


= s(B'  + ” T(w). 


Substituting  the  last  three  integrals  into  Eq.(83), 
we  have 

5 


T(w)  = s(P*  + iS-,p)w  _ sBw“^  _H(w)0'.''w)  + S~ A (85) 

p=o  ^ 


where 


H(w)  = (s  + tw^  + rw^')(sw^  _ 2tw^  _ hr)"^w^.  (86) 


Substituting  the  values  of  0(y/)  and  T(w)  from  Eqs.(82) 
and  (8^^)  into  Eq.(33),  we  have 


U = l^r(w)”'^  + t(w)"^  + swJ(sBw  + JZ 


) 


K(b. 

o 


n=l 


-1  . P 


+ 13^2)^  - sBw""^  + ^ A w 

p=o  ^ 


-(s  + tw^  + rw^)(sw^  _2tw^  _5r)“^(sBw^  _ ^na  w^"’^)^ 


(s  - 2tw“^  - '=;rw~^)dw, 


(87) 


where  the  real  nart  of  the  exnression  on  the  right 
is  to  be  taken.  The  value  of  the  constants  can  be 
obtained  from  Eqs.(84). 


CHAPTER  IV 


THE  STRESS  COMPONENTS 


In  the  preceeding  chapter  we  found  the  stress 
function  U in  terms  of  w.  With  the  aid  of  the  results 
obtained,  we  can  now  express  the  stress  components  in 
terms  of  w.  Prom  Eqs.(39)  and  (48),  we  have 


XX  + fy  = 4 Re(0»(z))  = 4 Re(|^-j^-^) 


(88) 


4 Re 


P*  (w)jgi*  (w) 
P» (w)P* (w) 


= 4M“  Re( P * (w) * (w) ) , 


where 


M s M(w)  =P*(w)P*(w),  (89) 

and  M is  real.  Prom  Eqs.(40),  (49),  and  (57),  we  have 


^(vy  - rx  + 2iQ)  » z0»(z)  + T*(z)  (90) 

_ P(w)P«  (w)jZl”(w)  . P(w)P**(w)gl»  (w)  ^ T'  (w) 

(P»(w))^  P*(w) 

= P(w)P»  (w)i?f"(w)  - P(w)gf*  (w)P"(w)  + P«  (w)P»  (w)T«  (w)  . 

(P’(i))-V 
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Let  us  denote  the  real  part  of  P*(w)j?f*(w)  by  J . 

0 

Prom  Eq.  (88) 

XX  + y>  = (91) 

Let  us  denote  the  real  part  of  (P'(w'l)^  by  J^,  and 
the  imaginary  part  by  If  we  denote  the  real  part 

of  P(w)pt(w)0"(w)  _ P(w)0.(w)P"(w)  + T*(w)p.(w)p.(w)  by 
J^,  and  the  imaginary  part  by  then  Eq.(90)  gives’ 

w M + ?i9  = 2ir’(j^  + iJjjKjj  + ijjj) 

(92) 

= 2M-’(JjJ2  . JjiJjj)  + 21M-’(JjJj2  + Vll’- 

Prom  Eqs.(91)  and  (92),  we  have 


" * - '■■•’(  V2  - JllJ22>’ 

(93) 

3^  = 2M-\  + 

(94) 

and 


Jy  = + Vii>- 


(95) 
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CASE  ONE. 

Let  r = 0,  and  let  t and  3 be  real. 

Setting  w = T?e^^,  we  have 

M = 8^  - 48tR"^Cos  3S  + 4t^R"^, 

P C 

= 8 B - 23tR”^(S^2^^^  ® B’Co8  B) 

- 8^R“^(S3^2^^”  - 4't^BR"'^, 

= 8^B*  + 23^tS^2^"^(S^”  ® 

- 28^tB*R”^(C03  B + C08  7B) 

+ 28BR”^(8^  + 2t^  _ 6t^R“^)Cos  2B 

- 4a^tP"^(B'  Cos  3B  + S3_2Sin  3B) 

+ sR~^(3a^  + 6t^  -2t^R“^  -2s^R^) (B’Cos  4B  -S^2^^” 

+ 4tBR"^(4t^  + 2s^  _ 3s^R^)Cos  5B 
+ 4st^R”^(B'C08  6B  + S^23^” 

= 3^  _ 6s^tR"^Cos  33  + 128t^R*So3  6B  - St^R’^Cos  9B, 

J22  = - 6s^tR-^Sin  33  + 12st^R~^Sin  6B  - 8t^R”'Sin  9B. 

r 
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J^l  = + 2a^tB’R"^(Sin  B + Sin  76) 

? C 

+ ?8 

- 28BR“^(s^  + 2t^  - 6t^R’^)Sin  26 

- 48^tR”^(S^2^o®  - B*Sin  36) 

-8R"^(3s^  ;|;  6t^  -2t^R*^  -23^R^)(B»Sin  46  + 

- 4tBR”'’(4t^  + 28^  - '3s^R^)Sin  5S 

+ 4st^R“^(S^2^°®  - B’Sin  66). 

These  values  of  F,  J^,  J^,  3 2*  ^22* 

substituted  into  Eqs.(93),  (94),  and  (95),  enable  us  to 
calculate  the  stresses  for  Case  I.  The  values  of  the 
stress  comoonenta  on  the  ed^^e  of  the  hole  of  Pig.l, 
Aonendix  I,  were  calculated  in  this  manner.  The  results 
are  given  in  Table  1,  Appendix  II . In  table  2,  Aopendix 
TI,  the  values  of  the  stresses  along  the  positive  x— axis 
exterior  to  f- e hole  are  given.  It  will  be  seen  that  as 
X increases,  the  stresses  rapidlv  approach  their  value 
at  infinity. 

The  stress  components  for  Case  IT  and  Case  III  can 
be  obtained  by  following  the  procedure  used  for  Case  I. 


CHAPTER  V 


SUMMARY 

In  this  dissertation  the  author  has  made  no  attempt 
at  exhausting  the  many  types  of  holes  obtainable  from 
Eq.  (4).  The  figures  given  in  Appendix  I represent  some 
of  the  simpler  holes  obtained  by  equating  certain  parame- 
ters to  zero.  If  the  reader  is  interested  in  the  stresses 
due  to  a certain  hole  included  in  Eq.(4),  then  any  un- 
necessary parameters  in  that  equation  should  be  set  equal 
to  zero.  In  such  cases  some  of  the  constants  given  in 
Eq.(84)  will  also  be  zero,  which  will  simplify  the  com- 
puting of  the  stresses  from  Eqs.(93)  to  (95). 

The  stresses  due  to  the  hole  given  in  Pig.  1,  can  be 

found  in  Appendix  II.  After  the  author  computed  these 

stresses,  T)r.  C.  B.  Smith  called  his  attention  to  a 
similar  hole 

X = Cos  2S  + 4 Cos  B, 

y = - Sin  2S  + 4 Sin  B, 

for  which  A.  E.  Green  has  computed  the  stresses.  His 
results  are  to  be  found  in  the  Proceedings  of  the  Royal 
Society  of  London,  Vol.  184,  pages  231-252.  The  method 
he  used  in  solving  his  problem  was  quite  different  from 
the  method  employed  here,  and  his  boundary  curve  was 
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X = naCos  S + b Cos  nB, 
y = - naSin  B + b Sin  nB. 

In  a paper  of  the  same  title  as  that  given  in 
reference  2 of  the  Bibliography,  Vartin  Greenspan 
determined  the  stresses  in  a large,  uniformly  loaded 
plate  weakened  by  a small  hole  of  an  approximately 
ovaloid  shape.  Greenspan  employed  a rather  laborious 
method  of  piecing  together  particular  solutions  of  the 
biharmonic  equation  for  the  stress  function  until  all 
the  boundary  conditions  could  be  satisfied.  This 
process  would  become  prohibitive  in  the  case  of  more 
complicated  boundary  conditions.  In  reference  2 of  the 
Bibliography,  Morkovin  applied  the  elegant  and  more 
general  method  of  N.  I.  MusheliSvili  for  solving  the 
same  problem.  Additional  information  on  MusheliSvili »s 
method  mav  be  found  in  references  7 and  8 of  the 
Bibliography. 


APPENDIX  I 


FIGURES  OP  THE  HOLES 


y = - Sin  20+5  Sin  0. 


35 


o o 


36 


y = - Sin  28+2  Sin  8 
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y 


= _ Sin  5S  + 20  Sin  8. 
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y 


sin  5S  + 10  Sin  3 
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y = Cos  'iS  - Sin  SS  + 10(Cos  S + Sin  S). 
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y = Cos  5B  - Sin  50  + 20(Cos  3 + Sin  S). 


PK 
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appendix  II 


'^able  1.  Values  of  the  stress  components 
on  the  edge  of  the  hole  given  In  figure  1 . 


s 

3^ 

0° 

0.000 

_ 1.000 

o.ooo 

10® 

.188 

_ .662 

- .038 

o 

o 

c\ 

- .132 

_ .206 

.284 

30° 

_ .214 

.16-=: 

.294 

c 

c 

- .303 

.220 

.488 

o 

o 

_ .410 

.252 

. 666 

60° 

- .495 

.286 

.8'=;7 

o 

o 

_ .673 

.325 

1.103 

o 

c 

CO 

_ .724 

.344 

1.478 

90° 

- .84  4 

.^36 

2.113 

100° 

- .709 

.020 

3.376 

110° 

.278 

.062 

4.440 

120° 

1.732 

.268 

3.732 

130° 

.822 

1.094 

.346 

140° 

.101 

_ .014 

- .146 

1*^0° 

- .104 

- .74'=; 

_ .041 

160° 

- .043 

- .Q34 

- .018 

170° 

_ .026 

- .086 

_ .0005 

180° 

0.000 

_ 1.^00 

0.000 
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Table  2.  Values  of  the  stress  coraponerts 

at  points  along  the  x-aris  exterior 
to  the  hole  f^iven  in  figure  1. 


X 

Q/S], 

6.  CO 

0.000 

_ 1.000 

6.33 

- .177 

- .325 

7.94 

.192 

.031 

10.23 

.543 

.035 

I'^.ll 

.770 

o 

o 

• 

20.06 

.871 

.019 

2'^.04 

.910 

.015 

30.03 

.925 

— 

.001 
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